Mott-Peierls Transition in the extended Peierls-Hubbard model 
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The one-dimensional extended Peierls-Hubbard model is studied at several band fillings using the 
density matrix renormalization group method. Results show that the ground state evolves from a 
Mott-Peierls insulator with a correlation gap at half-filling to a soliton lattice with a small band 
gap away from half-filling. It is also confirmed that the ground state of the Peierls-Hubbard model 
undergoes a transition to a metallic state at finite doping. These results show that electronic 
correlations effects should be taken into account in theoretical studies of doped polyacetylene. They 
also show that a Mott-Peierls theory could explain the insulator-metal transition observed in this 
material. 
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Since the discovery of the metallic phase of doped 
polyacetylene, this material has been extensively stud- 
ied , but the mechanism of the insulator-metal transi- 
tion observed upon doping is still poorly understood. R is 
known that both the Peierls instability and electronic cor- 
relations play a fundamental role in the formation and the 
properties of the insulating phase |^ and thus undoped 
polyacetylene is a Mott-Peierls insulator Therefore, 
ten years ago, Baeriswyl, Carmelo and Maki (s) proposed 
that the insulator-metal transition was also driven by 
the interplay of electron-electron and electron-phonon in- 
teractions. Within the restricted Hartree-Fock approxi- 
mation they have shown the possibility of such a Mott- 
Peierls insulator-metal transition in the Peierls-Hubbard 
model, which is the simplest model of polyacetylene in- 
cluding both interactions. Recently, several works us- 
ing sophisticated numerical many-body methods, such as 
the Gutzwiller variational wavefunction |^ and quantum 
Monte Carlo (QMC) simulations 0], have confirmed the 
occurrence of an insulator-metal transition in this model. 
On the other hand. Wen and Su, which have applied the 
the density matrix renormalization group (DMRG) tech- 
nique to this problem, have disputed the existence of this 
transition Hartree-Fock |^ and QMC Q simulations 
have also shown that the nearest-neighbor Coulomb re- 
pulsion opposes and can prevent the formation of a metal- 
lic state in the extended Peierls-Hubbard model. 

As an attempt to clarify this issue 1 have studied the 
properties of the extended Peierls-Hubbard model with 
parameters leading to a Mott-Peierls insulating ground 
state at half-filling. Accurate ground states and gaps 
are obtained for open chains up to 200 sites and different 
band fillings using the DMRG method and finite size 
effects have been carefully analyzed. Results show that 
the ground state of a Mott-Peierls insulator evolves to a 
soliton lattice upon doping. This soliton lattice is quali- 
tatively similar to the ground state predicted by simple 
electron-phonon models but both the gap and 

the amplitude of the lattice distortion decrease faster in 
the extended Peierls-Hubbard model than in these mod- 
els when the doping increases. An insulator-metal tran- 



sition occurs at a finite doping concentration in the ab- 
sence of nearest-neighbor electron-electron interaction, in 
agreement with previous studies of the Peierls-Hubbard 
model ||-|^. These results demonstrate that electronic 
correlations effects are important and should be taken 
into account in theoretical studies of doped polyacety- 
lene. They also confirm that a Mott-Peierls theory |^|^] 
could explain the insulator-metal transition observed in 
polyacetylene. 

The one-dimensional extended Peierls-Hubbard model 
is defined by the Hamiltonian 
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The operators c[„{cia) create (destroy) an electron of 
spin a at site i, nia — c^a^ia and ni = ni-\ -I- n£|. t 
is the resonance integral for an undistorted lattice and 
fixes the energy scale, A is the electron-phonon coupling 
constant, U and V are the on-site and nearest- neighbor 
Coulomb repulsion. As this model has an electron-hole 
symmetry, only hole doping is examined. The doping rate 
y is defined as the fraction of electrons removed from a 
neutral chain (which corresponds to a half-filled band). 
The usual dimerization order parameter describes the 
lattice degrees of freedom. A linear term with constant 
P is explicitly included in the lattice elastic energy in- 
stead of the constraint on the dimerization order param- 
eter used in previous works in order to reduce the 
average bond length variation. The value of P > is 
determined by the condition that the linear term in the 
elastic energy equals zero in the ground state configura- 
tion at half-filling. The lattice dynamics is completely 
neglected in this approach and the electron-phonon in- 
teraction is taken into account only through the coupling 
between electrons and a classical lattice relaxation. 
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To determine the ground state of the Hamihonian (|l|), 
one has to find both the lattice configuration {A^} and 
the electronic wavefunction which minimize the total en- 
ergy. Using a finite system DMRG algorithm one 
can compute the electronic ground state, its energy and 
the gradient of this energy (thanks to the Hellmann- 
Feynman theorem) for any given lattice configuration, 
and thus perform the minimization of the total energy 
with respect to lattice degrees of freedom {A^} [^,1). 
In principle, direct electronic excitations can also be ob- 
tained by calculating excited states of the electronic part 
of the Hamiltonian (|I]) for a fixed lattice configuration. 
Unfortunately, while the DMRG method gives excellent 
results for ground states, it is more difficult to obtain 
results for specific excited states. Therefore, I have only 
calculated charge gaps, which can easily be obtained from 
ground state energies for different band fillings p. The 
charge gap is believed to be equal to the lowest opti- 
cal absorption energy in the thermodynamic limit of the 
Peierls-Hubbard model {V = 0). In the extended Peierls- 
Hubbard model {V ^ 0) the relation between charge gap 
and optical gap is not known precisely but in this work I 
have assumed that both quantities are roughly equivalent 
in the thermodynamic limit. 

All calculations have been carried out for several chain 
lengths up to 200 sites and results have always been ex- 
trapolated to an infinite chain. Only open chains are 
considered because the DMRG method performs much 
better in this case than for periodic boundary conditions. 
Computations have been performed so that numerical er- 
rors on the ground state dimerization parameter Af are 
smaller than 10~^t. Numerical errors on gap values are 
estimated to be less than lO^^i at half filling and around 
10~^t away from half filling. All these estimations of 
the accuracy are based on an analysis of the behavior of 
DMRG results as a function of the number m of quan- 
tum states kept per block. The largest value of m used 
in this work ranges from 80 for short chains (50 sites) at 
half filling to 400 for long chains (200 sites) away from 
half filling. Truncation errors are typically between 10^^ 
and lO^*". I have also checked the accuracy of DMRG 
calculations against exact numerical results for long (up 
to 100 sites) non-interacting {U = V = 0) chains and 
against exact results for the one-dimensional Hubbard 
model {Ai = and U = 0) An excellent agreement 
has been found in both cases. 

My results at half filling are in good agreement with 
results obtained previously with DMRG and other 
many-body techniques like exact diagonalizations, quan- 
tum Monte Carlo simulations, and variational meth- 
ods They confirm that undoped polyacetylene is 
a Mott-Peierls insulator, which can be described with 
a reasonable accuracy by the extended Peierls-Hubbard 
models. I have determined appropriate parameters for 
polyacetylene by comparing model predictions to exper- 
imental values for the optical gap at half-filling, the 
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FIG. 1. Optical gap (in units of t) of the Peierls-Hubbard 
(PH) model, extended Peierls-Hubbard (EPH) model and 
SSH model as a function of doping. 

optical transition energies induced by photo-generated 
neutral and charged solitons and the neutral soliton 
spin density obtained from magnetic resonance exper- 
iments [|l|,^. This comparison shows that A = 0.1, 
U = AV = 2.5t and t = 2.7eV seem to be appropriate for 
polyacetylene in agreement with previous studies [^J|] . 

It is important to realize that in the Mott-Peierls 
regime the optical gap at half-filling is essentially a corre- 
lation gap as in the one-dimensional Hubbard model , 
although the electron-phonon coupling and the Peierls 
instability are responsible for features like the dimeriza- 
tion and the existence of solitons. For instance, DMRG 
calculations predict an optical gap Eg = 0.74t for the pa- 
rameters mentioned above. As the gap is only ^ O.OAt for 
U ^ V = and A = 0.1, electronic correlations account 
for at least 94% of the gap. Consequently, one expects 
this gap to be strongly reduced as soon as the system is 
doped because electronic correlations do not contribute 
to the formation of a gap away from half filling in the 
Hubbard model. A strong experimental evidence for this 
reduction is the difference between the gap at half filling 
(1.8eT^) and the energy of the optical transition induced 
by photo-generated charged solitons {0A5eV) which cor- 
responds to the gap of a lightly doped chain in our sim- 
plified model. 

I have investigated the extended Peierls-Hubbard 
model at several dopant concentrations up to y = 16% 
in the Mott-Peierls regime. In this regime the system 
evolves upon doping from the Mott-Peierls insulating 
state to a soliton lattice with a small gap. The evolu- 
tion of the optical gap upon doping is shown in Fig. gfor 
the polyacetylene parameters mentioned previously. As 
expected, the gap is strongly reduced to ^ 0.17t as soon 
as the system is doped. The amplitude of the lattice dis- 
tortion and the gap decrease with increasing doping but 
no transition to a metallic state is found up to the highest 
doping studied in this work {y = 16%). It is possible that 
a transition occurs at a higher doping but this would not 
be relevant for the transition observed in polyacetylene 
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FIG. 2. Charge gap (in units of t) of the Peierls-Hubbard 
model as a function of the inverse system size for two different 
doping levels. Lines are linear extrapolations. 

around y — 6%. Quantum Monte Carlo simulations |^] 
have also shown that the lattice distortion survives at 
high doping for V = U/2. For y > 4%, the amplitude of 
the lattice distortion Ae corresponds exactly to the value 
of the gap if both quantities are extrapolated to an in- 
finite chain. Therefore, away from half filling the gap is 
a band gap generated by the lattice modulation, though 
electronic correlations contribute indirectly to its forma- 
tion because they increase the amplitude of the lattice 
distortion ||^. 

The soliton lattice found in the doped extended Peierls- 
Hubbard model is qualitatively similar to the soliton lat- 
tice predicted by simple electron-phonon models [p^ . 
However, it is important to realize that the evolution 
of the gap and lattice distortion amplitude upon dop- 
ing is different from the predictions of the Su-Schrieffer- 
Heeger (SSH) model (shown in Fig. || for A = 0.2). To 
reproduce the results obtained in the extended Peierls- 
Hubbard model with the SSH model one would need 
to use an effective electron-phonon coupling which de- 
creases with increasing doping. Moreover, such a doping- 
dependent parameter A(t/) should change abruptly at 
half-filling to reproduce the sudden disappearance of the 
correlation gap. Finally, one notes that the energy scales 
involved in both approaches differ by an order of mag- 
nitude. For instance, my calculations indicate a gap of 
about 0.07 eV at y = 8% while the SSH model pre- 
dicts OAeV These results demonstrate that elec- 
tronic correlations effects in such one-dimensional sys- 
tems are not reproduced by simple single-electron mod- 
els with an effective electron-phonon coupling, contrary 
to a basic assumption of the SSH theory of conducting 
polymers |p],pT|. Thus, the electron-electron interaction 
and electronic correlations effects should be taken into 
account in theoretical studies of doped polyacetylene. 

In Fig. 1^ one can also see that the gap of the Peierls- 
Hubbard model with A = 0.1 and U = 2.5t vanishes at a 
critical doping between 8 and 12%. These results confirm 
the existence of a transition to a metallic state at a finite 
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FIG. 3. Lattice dimerization 
Peierls-Hubbard model at y = 12%. 



parameter Am of the 



doping They also show that the electron-electron 

interaction can either support or oppose the Peierls in- 
stability away of half-filling depending on the parame- 
ters used, as predicted by the restricted Hartree-Fock 
approximation |^,^. It should be noted that my numer- 
ical results agree quantitatively with those presented in 
Ref. 1^ but additional calculations and an analysis of 
finite size and chain edge effects lead to a different con- 
clusion. Fig. H shows the charge gap as a function of 
the system size for y = 8% and 12%. The value of the 
gap extrapolated to an infinite chain is clearly finite for 
y = 8% but vanishes for y = 12% within numerical er- 
rors (~ 10~^t). The lattice dimerization parameter 
is shown in Fig. ^ for a 200-site chain at y = 12%. The 
shape of looks similar for y = 8% except that the 
amplitude is smaller at higher doping. However, in the 
insulating phase {y < 8%), the amphtude of the distor- 
tion in the middle of the chain tends to a finite value as 
the chain length increases. This confirms that this lattice 
modulation is a genuine Peierls distortion. On the other 
hand, in the metallic regime {y > 12%) the amplitude de- 
creases as a power-law with an exponent —0.66 as the sys- 
tem size increases. This behavior demonstrates that the 
lattice distortion is a chain edge effect with a very slow 
asymptotic decay. The amplitude of these Friedel oscilla- 
tions is known to decrease asymptotically as a power-law 
with an exponent -1 in a one-dimensional Fermi liquid, 
but density fluctuations are strongly affected by electron- 
electron interaction and the exponent can be as small as 
— i in a Luttinger fiquid |T5| . 

In conclusion I have investigated the ground state of 
the extended Peierls-Hubbard model in the Mott-Peierls 
regime, which is appropriate for polyacetylene. Results 
show that the ground state evolves from a Mott-Peierls 
insulator with a correlation gap at half-filling to a soliton 
lattice with a small band gap away from half-filling. It 
is also confirmed that a insulator-metal transition occurs 
in the Peierls-Hubbard model at a doping concentration 
between 8 and 12%. These results clearly show that elec- 
tronic correlations effects are important and should be 
taken into account in theoretical studies of doped poly- 
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acetylene. They also suggest that the primary mecha- 
nism of the insulator-metal transition in polyacetylene 
is the interplay between electron-phonon and electron- 
electron interactions, which induces a transition upon 
doping from an insulating state with a gap of ~ l.SeV 
to a state with a gap which is two orders of magnitude 
smaller. Obviously, this theoretical investigation of the 
properties of an ideal, infinite and isolated chain cannot 
describe the properties of actual physical systems. Un- 
derstanding the insulator-metal transition of polyacety- 
lene will require the study of more realistic model in- 
cluding lattice dynamics, interchain couplings, interac- 
tion with dopant ions and disorder |p^ . In the future I 
plan to study such models with DMRG, particularly the 
effects of quantum lattice dynamics. 
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